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ABSTRACT 

We investigate the latitude at which type I X-ray bursts are ignited on rapidly rotating accreting neutron 
stars. We find that, for a wide range of accretion rates M, ignition occurs preferentially at the equator, in accord 
with the work of Spitkovsky et al. However, for a range of M below the critical M above which bursts cease, 
ignition occurs preferentially at higher latitudes. The range of M over which nonequatorial ignition occurs 
is an increasing function of the neutron star spin frequency. These findings have significant implications for 
thermonuclear flame propagation, and they may explain why oscillations during the burst rise are detected pre- 
dominantly when the accretion rate is high. They also support the suggestion of Bhattacharyya & Strohmayer 
that non-photospheric radius expansion double-peaked bursts and the unusual harmonic content of oscillations 
during the rise of some bursts result from ignition at or near a rotational pole. 
Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries — X-rays: bursts 
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1. INTRODUCTION 

X-ray bursts are thermonuclear explosions 
on the surfaces of accreting neutron stars in 
X-ray binarie s (LM X Bs) jBa bushkina et all 
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they are triggered b y unstable hydrogen or helium burning 
(for r eviews, see Cumming 2004; Strohmaver & Bildsten 
2006). The physics of type I X-ray bursts is generally well 
understood, and detailed time-dependent one-dimen sional 
models (e.g., [Woosley et al. 2004; iFisker et"aTJ 120061) have 
been rather successful at reproducing the gross characteristics 
of bursts, such as their fast rise times of ~ 1 s, durations of 
~ 10-100 s, and recurrence times of a few hours to days. 

Nearly all type I X-ray burst models are one-dimensional 
and thus implicitly assume that matter accretes spherically 
onto a nonrotating neutron star. Therefore, ignition occurs 
simultaneously over the entire stellar surface. However, these 
assumptions are clearly inapplicable to the vast majority of 
accreting neutron stars. The Rossi X-Ray Timing Explorer 
has detected highly sinusoidal oscillations with frequencies of 
45-1 122 Hz during bursts in 17 L MXBs (Strohmayer & Bild- 
sten 2006 and references therein; rBhattacharvva et aLll2006l : 
iBhattacharvval 120061; iKaaret et al.1 120061) . It is thought that 
the burst oscillation frequency corresponds to the neu tron star 
spin frequency (e.g. Strohmaver & Mar kwardtll 1999b . which 
implies that many neutron s tars that exhibit bursts a re in 
fact rapidly rotating. Indeed, Mu no et all (1200 lL l2004h have 
shown that the properties of bursts depend sensitively on the 
neutron star rotation rate, and so rotation must be consid- 
ered in theoretical burst models. Furthermore, since the time 
required to accumulate a sufficient amount of fuel to trig- 
ger a burst is much greater than the duration of the result- 
ing burst, it is highly unlikely that all of the accreted fuel 
over the en t ire neutron s tar surface will ignite simultaneously 
dJosslll978l ; ISharal[l982l) . Thus, contrary to the assumptions 
of most theoretical models, ignition almost surely occurs at a 
point, and the resulting thermonuclear flame subsequently en- 
gulfs the whole stellar surface in ~ 1 s, the burst rise time 
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dFrvxell & Wooslevl 119821; iBildstenl [l99l ISpitkovskv etail 
120021 hereafter SLU02). Time-resolved spectroscopy and ob- 
servations of large amplitude oscillations during the rise of 
some type I X-ray bursts from rapidly rotating neutro n stars 
support this conclusion (Strohmave r et alJll997l lT998). 

SLU02 investigated the dynamics of localized ignition on 
a rapidly rotating neutron star. They found that the latitude 
at which localized ignition occurs is of great importance in 
determining the burning front propagation speed. It can also 
affect the stability a nd nature of zon al flows during the de- 
cay phase of bursts (Cumming 20 051), the harmonic conten t 
of burst rise oscillations (Bhattacha ryya & Strohmaverl 2005), 
and the light curves of non-photospheric radius expansion 
(PRE) bursts dBhattacharvva & Strohmaverl2006allb1) . SLU02 
assert that, due to the reduction of the effective gravita- 
tional acceleration, ignition is likely to occur at the equator 
of a rapidly rotating neu tron star. However, recent observa- 
tions s uggest otherwise. Bhattacha ryya & Strohmaverl d2005l 
2006a b) argue that the harmonic content of burst rise oscilla- 
tions and non-PRE double-peaked bursts from the LMXB 4U 
1636-536, which harbors a neutron star with spin frequency 
v = 581 Hz, require ignition to oftentimes occur at or near the 
rotational pole, at least in this source. 

In this Letter, we evaluate theoretically the latitude at which 
type I X-ray bursts on rapidly rotating neutron stars are most 
probably ignited. In $2] we outline the physics that governs 
the latitude of type I X-ray burst ignition and determine the 
ignition latitude as a function of the global accretion rate. We 
discuss our results and conclude in $3] 

2. IGNITION LATITUDE 

In this section, we first outline the basic physics that deter- 
mines the latitude of type I X-ray burst ignition on a rapidly 
rotating neutron star using a simpl e model. We then use the 
global linear stability analysis of Cooper & Nara yanl d2005l) 
to carry out more detailed calculations. 

2.1. Basic Physics 

Using the one-zone burst ignition model of Bildsten ( 1998), 
we expand upon the work of SLU02 to determine the most 
probable latitude of burst ignition on the surface of a rapidly 
rotating neutron star. The reader is encouraged to refer to 
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these works for further details. We neglect general relativis- 
tic corrections throughout for clarity. Following SLU02, we 
assume that at all times prior to ignition the accreted plasma 
is in hydro static equilibrium and is at rest in the corotating 
frame (e.g. Inogamov & Sunyaev 1999). We presume that the 
accreted matter spreads in such a way as to minimize the grav- 
itational potential energy of the accreted layer. The base of 
the accreted layer is therefore an equipotential surface, and so 
the pressure at the base of the acc reted layer is the same at 
every latitude (e.g., Clayton 1983). Hydrostatic equilibrium 
thus implies that P = £i ayer (A),g e ff(A) is independent of lati- 
tude, where P is the pressure at the base of the accreted layer, 
Siay er is the column depth, g e g is the effective gravitational ac- 
celeration, and A = 7r/2-6>is the latitude. Thus, £i ay er oc g~^, 
and so 

Socgeff, (1) 

where £ = £(A) is the local ma ss accreti o n rate per unit area. 
According to equation (20) of Bildsten] (1 1998b . the column 
depth at which helium ignites E; gn varies as 
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Combining equations ([TJ and (f2]i gives the ignition time as a 
function of latitude 



£ign/£0Cgef f 5 . 
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Thus fi gn , the time it takes to accrete a critical amount of fuel 
to trigger a type I X-ray burst, is an increasing function of 
g e ff. Due to centrifugal acceleration, g e s is lowest at the equa- 
tor of a rapidly rotating neutron star, which means that ?i gn (A) 
is a minimum at the equator. Therefore, ignition will occur 
preferentially at the equator (SLU02). 

However, nuclear burning on the surface of an accreting 
neutron star does n ot always trigger a burst. Both theo- 
retical models (e.g.. iFuiimoto et all 119811; iPaczvnskH 1 1 983b 
Bildste n 1998) and o bservations (e.g. ICornelisse et alj|2003t 
Galloway et al.i r2006) suggest that there is a local critical ac- 
cretion rate E^t above which nuclear burning is stable, and 
thus bursts do not occur. By equation ((TJ, £ cr j t is related to 
the global critical accretion rate M such that 
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where M cl -i t = M cl -i t ( A) is defined to be the critical global accre- 
tion rate above which nuclear burning is stab le at latitude A. 
According to equation (24) of Bildsten ( 1998), 
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This implies that, since g e a is lowest at the equator and in- 
creases towards the rotational pole, there exists a range of 
accretion rates for which nuclear burning is stable near the 
equator and unstable near the pole. Consequently, for a given 
Ao and corresponding critical accretion rate M cr ; t (Ao), nuclear 
burning is stable for all A < Aq and unstable for all A > Ao. 
Clearly, at these accretion rates, bursts will be triggered off 
the equator. For M > M cl -i t (7r/2), no bursts are triggered at 
any latitude. Taking the functional derivative of equation © 
gives 
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FIG. 1 . — Plot shows Aj gll , the latitude at which type I X-ray bursts ignite on 
a rapidly rotating neutron star, derived from a one-zone burst model for two 
stellar rotation frequencies v. M = 1 .4Mq and R = 10 km. For a wide range of 
global accretion rates M, ignition occurs at the equator (A = 0). However, for 
a range of M near the critical rate above which bursts cease, ignition occurs 
at higher latitudes. The range of M over which nonequatorial ignition occurs 
is oc v 2 . 

where v, M, and R are the spin frequency, mass, and radius 
of the neutron star, respectively. Thus, for a neutron star with 
spin frequency v = 600 Hz, type I X-ray bursts will ignite at 
the equator forM/M cr i t (7r/2) < 0.89 and off of the equator for 
0.89 <M/M crit (7r/2) < 1. Within the latter range of accretion 
rates, the latitude at which ignition occurs is given by 
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where vk_ is the Keplerian frequency. Figure Q] shows a plot 
of Aign as a function of M. We see that type I X-ray burst ig- 
nition on a rapidly rotating neutron star occurs preferentially 
near the equator for a wide range of M. However, for a small 
range of M near M c ,i t , ignition occurs at higher latitudes be- 
cause nuclear burning becomes stable near the equator. As 
equation (|7]) illustrates, the range of accretion rates over which 
nonequatorial ignition occurs increases dramatically with v. 

2.2. Global Linear Stability Analysis 

We now use t h e glo bal linear stability analysis of 
ICooper & Naravanl d2005l) . which is an expanded a nd im- 
proved version of the model of Narayan & Heyl (2003), to de- 
termine the latitude of type I X-ray burst ignition on rapidly 
rotating neutron stars. We assume that matter accretes at a 
global rate M onto a rapidly rotating neutron star of mass 
M = 1.4M0, radius R = 10 km, and spin frequency v = 650 
Hz, and it spreads over the stellar surface in the same manner 
as described in 32.11 We set the composition of the accreted 
matter to be that of the Sun, such that at the neutron star sur- 
face the hydrogen mass fraction X = 0.7, helium mass frac- 
tion Y = 0.28, CNO mass fraction Zcno = 0.016, and heavy 
element fraction Z = 0.004, where Z refers to all metals other 
than CNO. 

We plot the ignition time fj gn as a function of M for A = 0, 
7r/4, and 7r/2 in Figure |U For M < 0.16M Edd , f ign (A) is a 
minimum at the equator, and so bursts ignite preferentially 
near the equator, in agreement with both SLU02 and equa- 
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FIG. 2. — Plot shows fi sn (A), the type I X-ray burst ignition timescale, 
at three different latitudes A for a neutron star with spin frequency v = 650 
Hz, M = 1.4M©, and R = 10 km. For accretion rates M < 0.16MEdd, ?ign 
is lowest at the equator, which means that bursts ignite at the equator. For 
0.16 < M/MEdd ^ 0.19, the burst ignition latitude increases with M. For 
0.19 < M/MEdd 0-3, although nuclear burning is unstable at all A, ignition 
occurs preferentially near the pole where t\ gn is lowest. For 0.3 < M/MEdd ^ 
0.34, nuclear burning becomes stable at low latitudes, and ignition still occurs 
preferentially near the pole. 

tion (0). Here, Msdd = 8.3 x 10 17 gs" 1 denotes the mass ac- 
cretion rate at which the accretion luminosity is equal to the 
Eddington limit. For 0.16 < M/M Edd < 0.19, A ign gradually 
increases with increasing M, again in accord with the results 
of j j2.1| However, Figure [2] shows a new regime of bursts 
for M > 0.19MEdd> the regi me of "delayed mixed b ursts," 
which was first identified by Nar avan & Hevll (|2003). De- 
layed mixed bursts are mixed hydrogen and helium bursts trig- 
gered by unstable helium burning, and they are preceded by 
a long period of stable nuclear burning. Using their global 
linear stability analysis, Narayan & Heyl found that, within 
this range of accretion rates, (1) the ignition timescale is an 
increasing fu nction of M (which is consistent with obser- 
vations, e.g., [van Paradiis et al.lll979lll988t [Cornelisse et al.l 
l200l iRerrrillard et al.1 l2006t iGallowav et all 120061) . and (21 
the critical global accretion rate above which delayed mixed 
bursts occur is an increasing function of the gravitational ac- 
celeration geff- Thus, for the range 0.19 < M/MEdd < 0.3, 
although nuclear burning is unstable at all latitudes, ignition 
occurs preferentially near the pole where f; gn is lowest. For 

0.3 < M/MEdd ^ 0.34, nuclear burning becomes completely 
stable at low latitudes, and ignition again occurs preferentially 
near the pole. For M/MEdd ^ 0.34, nuclear burning is stable 
at all latitudes, and no bursts occur. 

3. DISCUSSION AND CONCLUSIONS 

The results of both the simple one-zone model of lBildstenl 
(119981) and the more detaile d global linear stability analysis 
of ICooper & Naravanl (120051) suggest that bursts ignite on the 
equator at low accretion rates and off of the equator at higher 
accretion rates. The global linear stability analysis predicts 
that nonequatorial ignition occurs over a much larger range 
of M - nearly 50% for v = 650 Hz versus only about 10% 
in the one-zone model - and that ignition should occur near 
the pole for most of this range. The differences in the re- 
sults of the two models are due to the delayed mixed burst 
regime, which is predicted in the global linear stability anal- 
ysis but not in the one-zone model (for an explanation of 



delayed mixed bursts, see ICooper & Na ravan 20061) . In ef- 
fect, the delayed mixed burst regime significantly extends the 
range of M over which nonequatorial ignition occurs. The 
global linear stability an alysis and the two-zone model of 
Cooper & Naravan (2006) agree better with observations than 
all other current burst models, but this agreement h olds only 
in a time-averaged sense (e.g. Gal loway et al.ll2006h . By this 
we mean that, while observations imply that the mean ignition 
time (fj gn ) is an increasing function of M for M > 0.15M£dd, 
which is in accord with the results shown in Figure [2] the 
ignition time ?i gn measured between pairs of bursts observed 
in nature exhibits significant deviations relative to the mean. 
Our models cannot account for this behavior. Since the pre- 
dicted burst ignition latitude is model-dependent, and no cur- 
rent model can successfully reproduce the chaotic behavior of 
bursts, we have little confidence in the accuracy of our pre- 
dicted ignition latitudes A; gn . However, all theoretical models 
would predict that bursts ignite preferentially near the equator 
at low values of M and off of the equator at higher M. This 
result is very robust. 

Recent observations have generated a renewed interest in 
the latitude of type I X-ray burst ignition. First, non-PRE 
doubl e peaked bursts have been observed in a few L MXBs 
(e.g. [S ztaino et aTlll985b iPenninx et al]|1989t; [Kuulkers et alj 
120021) . IBhattacharyya & Strohmayerl (l2006albf) argue that ig- 
nition near a rotational pole can explain the non-PRE dou- 
ble peaked bursts observed in 4U 1636-536, which contains 
a neutron star with v = 581 Hz. Their simple model of 
nonequatorial ignition and subsequent thermonuclear flame 
propagation and temporary stalling near the equator quali- 
tatively reproduces both the light curves and spectral pro- 
files of such bursts. Second, although burst oscillations de- 
tected in both the rise and decay are usually quite sinusoidal, 
IBhattacharyya & Strohmaveii d2005l) report evidence for sub- 
stantial harmonic content in the oscillations during the rise of 
a burst from4U 1636-536. They again suggest that nonequa- 
torial ignition and subsequent flame propagation can explain 
the observations. These authors acknowledged that nonequa- 
torial ignition is at odds with the work of SLU02, who argued 
for equatorial ignition at all M. Our result that nonequatorial 
ignition is more likely when the accretion rate is high natu- 
rally reconciles this discrepancy. 

That nonequatorial ignition is preferred at higher M may 
explain why oscillations in the rising phase of some bursts 
on rapidly rotating neutron stars occur predominantly on the 
banana branch, i.e. when the inferred accretion rate is high 
(IStrohmayer & Bildstenl20 06). This can be understood as fol- 
lows. SLU02 find that burning fronts propagate much faster 
near the equator than near the poles. If M is low and ignition 
occurs near the equator, the flame will spread in longitude on a 
timescale that is much shorter than the burst rise timescale and 
quickly produce an axisymmetric belt. An observer would 
detect no oscillations during the rising phase because there is 
no longitudinal asymmetry. On the other hand, if M is high 
and ignition occurs close to, but not directly at, a rotational 
pole, the flame will spread in the longitudinal direction on a 
timescale that is comparable to the burst rise timescale. This 
slow longitudinal spread would create a non-axisymmetric hot 
spot and hence oscillations. We note that oscillations in the 
decay phase of some bursts are again observed only when a 
system is on the banana branch. However, burst oscillations in 
the decay phase are not well understood theoretically, and we 
are hesitant to suggest a connection between these oscillations 
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and nonequatorial ignition. 

We stress that our results apply only to rapidly rotating, 
weakly-magnetic neutron stars. Rapid rotation induces a 
strong latitudinal dependence on the effective gravitational ac- 
celeration which restricts burst ignition to a preferred latitude. 
For a nonrotating (or slowly rotating) neutron star, all physi- 
cal quantities are independent of A (or nearly so), and ignition 
could occur at any latitude. Thus, oscillations during the burst 
rise could occur at any M for slowly rotating neutron stars. 

Strong magnetic fields may channel and confine accreted 
matter to some region of the neutron star surface. In this case, 
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